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Abstract. We calculate the thermochemical equilibrium 
of the diffuse interstellar medium, including ionization by 
a photon flux Fi, from neutrino decay. The main heating 
mechanism considered is photoelectrons from grains and 
PAHs. For the studied range of F,^ values, there always 
exists two regions of stability (a warm and a cold phase) 
that can coexist in equilibrium if the thermal interstellar 
pressure is between a maximum value {Pmax) and a mini- 



raax } 

mum value (P™„). High F^ values (- lO'^-lO^ cm'^ s'^) 
can be consistent with observed interstellar pressures only 
if more efficient sources are heating the gas. It is shown 
that a neutrino flux increase (due, for example, to an in- 
crease in the supernova explosion rate) may stimulate the 
condensation of cold gas by decreasing Pmax below the 
interstellar pressure value. 
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1. Introduction 

The diffuse interstellar medium (ISM) is observed to be 
inhomogeneous with cold (T ^ 10^ K) clouds embedded 
in a warmer (T ~ 10** K) intercloud gas (see, for example, 
the review of Kulkarni & Heiles 1988). The theoretical 
explanation for this structure was provided by Field et 
al. (1969) who showed that two thermally stable phases 
can coexist in pressure equilibrium over a limited range 
of pressures, close to those observed in the ISM. Since the 
two-phase model of Field et al. (1969), the thermal and 
ionization equilibrium of the ISM and its stability have 
been studied by many authors including different heat- 
ing and ionizing processes (Black 1987; Kulkarni & Heiles 
1988). The response of the ISM to variations in physical 
processes and parameters are of interest both to a better 
understanding of the ISM behavior and because of its cen- 
tral role in star formation and galaxy evolution models. Ef- 
fects on the ISM equilibrium of variations in, for example. 
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X-ray and far UV radiation fields, cosmic ray ionization 
and metal abundance have been studied by several authors 
(ShuU & Woods 1985; Parravano 1987; Wolfire et al. 1995; 
Parravano & Pech 1997). In this work, we are interested 
in a particular process: the flux of ionizing photons com- 
ing from the radiative decay of neutrinos (Sciama 1990). 
Notwithstanding this is a speculative theory, Sciama has 
argued in a set of papers (Sciama 1993, 1995, 1997a, 
1997b, 1998) that it can explain the widespread ionization 
far from the galactic disk and many other observational 
results. The simplification in Sciama's work is to assume 
a temperature of ^ 10* K without explicitly solving the 
thermal equilibrium. The former was made by Dettmar 
& Schulz (1992) who showed that heat input associated 
with neutrino decay is too small to account for the ob- 
served ISM temperature. However, their conclusion that 
neutrino decay cannot be a dominant source of ionization 
could be mistaken, as was pointed out by Sciama (1993), 
because they neglected the existence of other known heat- 
ing mechanisms. Our goal here is to analize, in a more 
complete and self-consistent way, the effect of ionization 
due to neutrino decay on the thermochemical equilibrium 
of the ISM. 

In Sect. 2 we discuss the physical processes included 
in this work and provide the basic equations. Sect. 3 is 
dedicated to a discussion of the results, and the main con- 
clusions are summarized in Sect. 4. 



2. Basic equations 

In order to analize the effect of neutrino decay photons 
on the thermochemical equilibrium, a simple model for 
the ISM is used. The included cooling mechanisms are: a) 
cooling by collisions of electrons with 0+, Si"*", Fe"*", 0+, 
S+ and N (Ag); b) cooling by collisions of neutral hydro- 
gen with Si+, Fe+ and C+ (A//); and c) cooling due to 
Ly-a excitation by electrons (A^j,). All the cooling rates 
and the relative abundances were taken from Dalgarno & 
McCray (1972). We consider the following heating mech- 
anisms: 

a) Interaction of cosmic ray with hydrogen atoms and elec- 
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trons: 

Fcr = Ccr [5 X 10-12(1 + $)n(l - x) + 5.1 X 10-i"ne] 

ergs cm-^ s""^ , (1) 

where n = n(HI) + n(HII) is the total number density 
of hydrogen, Ue is the number density of electrons and 
X = n(HII)/n is the ionization degree of hydrogen. The 
number of secondary ionizations ($) was taken from Dal- 
garno & McCray (1972), and the primary ionization rate 
is assumed to be C,cr = IQ--^'' (Spitzer 1978; Black et 
al. 1990; Webber 1998). 

b) Heating by H2 formation on dust grains (Spitzer 1978): 

= 4.4 X lQ-'^^n^{l - x) ergs cm"^ s^^ . (2) 

c) Photoelectric heating from small grains and PAHs 
(Bakes & Tielens 1994): 

Fpe = IQ-^^eGoTi ergs cm"^ s"^ , (3) 
where the heating efficiency (e) is given by 
_ 4.87 X 10"- 

[1 + 4 X 10-3 (GoTV2/ne)°-^^l 

3.65 X 10-2 (r/104)°-^ 
[1 + 2 X 10-4 (C?„Ti/2/n,)] ' ^"^^ 

and Go is the far UV field normalized to its solar neigh- 
borhood value. 

We only consider ionization/recombination for hydro- 
gen. The recombination rate is given by 

X+ = xnea{T) s'^ , (5) 

where a{T) is the recombination coefficient to all states 
except the ground one, and it was taken from Spitzer 
(1978). The rate of ionization by cosmic rays is given by 

X~ = Ccril + m - X) . (6) 

We also use the simple analytic fits provided by Wolfire et 
al. (1995) to estimate the ionization {X^j^) and heating 
(Txb) due to the soft X-ray background as functions of 
the column density (Ny^) and the electron fraction in^jn). 
In this work we adopt = lO^^cm-^. 

In addition to the above sources of ionization, we also 
consider the photons produced by neutrino decay. The ion- 
ization due to this mechanism can be written in the form 
(Sciama 1990): 

X- = F,a(l - x) , (7) 

where a = 6.3 x 10- cm^ is the absorption cross sec- 
tion of hydrogen and F^, is the flux of hydrogen-ionizing 
photons produced by neutrino decay. In this work F^, is a 
free parameter, although Sciama (1997a) estimated that a 
value of ~ 3 X 10'' photons cm-^ is necessary to pro- 
duce an electron density rig ~ 0.05 cm-^ in the intercloud 
medium. The most recent (but still uncertain) estimation 
was F^ ^ 10^ photons cm-^ (Sciama 1998). 



3. Results and discussion 

The thermochemical equilibrium is calculated by solving 
simultaneously the equations A = F and X+ — X~, 
where A = Ae -I- Ah + A^j, is the total cooling rate, 
r = Per -I- Th + Ppe + Tjcfl is the total heating rate, 
and X~ = X~j. + X^j^ + X~ is the total ionization rate. 
Fig. la shows the equilibrium pressure-density relations 
for Go = 1 and for three different values of F^, (0, 10^ 
and 10'' cm-2 s-^). The corresponding electron fractions 
(ne/n) are showed in Fig. lb, where it can be seen that, as 
expected, Ue/n increases as F^ increases. Most of the ion- 
ization for the cases F^, > 10 cm-^ s-^ is due to photons 
coming from neutrino decay, and thus the neutrino decay 
is a very efficient ionization mechanism. An increase in 
Fi, (and the consequent increase in the electron density) 
enhances the cooling by electron collisions (Ag). Addition- 
ally, the dominant heating mechanism is always photoelec- 
trons from grains and PAHs (Fpe ) , which almost it is not 
affected by the flux F^. Consequently, for a given den- 
sity, when Fi, increases the thermal equilibrium is reached 
at lower temperatures in the regions where Ag dominates 
the cooling (high densities), and the equilibrium curve is 
shifted down (see Fig. la). In contrast, at low densities 
(n ~ 10-2 cm-3), the dominant cooling mechanism is 
A^y, which decreases when F^, increases and, in this case, 
the equilibrium is reached at higher temperatures. 

Fig. la also shows that two regions of thermal stabil- 
ity, i.e., where the slope is positive (Field 1965), always 
exist and two phases can coexist in pressure equilibrium 
if the interstellar pressure P = p/k is between a mini- 
mum {Pmin) and a maximum (Pmax) value. For the case 
F„ = we obtain Pmax — HOO K cm-^ and Pmin — 490 
K cm-3; and if we assume an equilibrium pressure of 
P ~ 10^ K cm-3, then there can be gas with T ~ 8800 
K, n ~ 0.1 cm-3 and rie/n 4.7 x 10-^, and gas with 
T ~ 100 K, n ~ 10 cm-^ and ^ 1-7 x 10-^ co- 

existing in equilibrium. These results agree roughly with 
observational estimations of the warm and cold neutral 
phases in the local ISM (Kulkarni & Heiles 1988). Fur- 
thermore, it can be seen in Fig. la that Pmax, the max- 
imum pressure value over which only the cold phase can 
exist, decreases as F^, increases. For F^ = 10^ cm-^ 
and for F^ = 10'' cm-^ s-^ we obtain Pmax ~ 640 K 
cm-3 and Pmax ~ 500 K cm-^, respectively; but observa- 
tions indicate that the pressure in most of the regions of 
the Galactic plane is ^ 10"^ K cm"'^ (Jenkins et al. 1983). 
Therefore, there seems to be an inconsistency between the 
observed pressure in a multi-phase medium and high 
values. The basic reason for this behavior (that high F^, 
values imply too low Pmax values) is that neutrino decay is 
a poor heating agent, while other processes can ionize and 
heat the ISM. For instance, when the cosmic ray ioniza- 
tion rate (Ccr) is changed (keeping F^ = fixed), we find 
that Pmax decreases as Ccr increases, but for Ccr 10-"'^ 
the heating by cosmic rays (Per) becomes more im- 
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portant that heating by photoelectrons from grains and 
PAHs (Fpe), and then Pmax begins to increase as Ccr is 
increased. The minimum Pmax value is ^ 900 K cm^'^. 

However, more efficient heating mechanisms that those 
considered here can rise up the equiUbrium curve increas- 
ing Pmax- In order to illustrate this effect, we have plotted 
in Fig. 2 Pmax as a function of F^, for three different values 
of Go (1, 10 and 20) and, in consequence, three different 
values of Tpe (believed to be an important heating mech- 
anism in the ISM). We can see that as F^, is increased 
Pmax decreases until a minimum value ('^ 480 K cm~'^ 
for Go = 1, - 590 K cm"^ ^j. ^ iq j^^^j _ ggO K 
cm~^ for Go = 20) and after that remains constant. This 
occurs when Ue/n — > 1 in the warm gas and, therefore, ad- 
ditional increases in F,y do not produce additional changes 
in this phase. Fig. 2 shows that if an ISM pressure of ~ 10"^ 
K cm^"^ is assumed, a two-phase medium is possible for 
Go = 1 only if F^ <, 10 cm^^ g-i, and for Go = 20 only if 
Fij ^ 200 cm~^ s~^. We conclude that high fluxes of neu- 
trino decay photons 10'^ cm~^ s~^) can be consistent 
with a two-phase medium only if more efficient heating 
sources are acting on the gas. 

Sciama (1997a) estimated that F^ ^ 3 x 10^ cm~^ s~'^ 
is necessary to produce Ue ~ 0.05 cm~'^ at T ~ 6000 K. 
Fig. 3 shows Tie as function of F^, for T = 6000 K and for 
the same three values of Go given in Fig. 2. The desired 
electron density at this temperature can be reached at 
high values only if Go ^ 20 and, again, more efficient 
heating sources seem to be necessary. 

An interesting consequence has to be noted: neutrino 
beams from neighboring regions may induce the conden- 
sation of cold clouds stimulating the formation of stars. 
The importance of star formation triggered by previously 
formed stars has been recognized by many authors (see 
the review of Elmegreen 1992). Triggering mechanisms are 
usually related to compression of the ISM by shock waves 
from close supernovas, because the transition warm gas 
cold gas is promoted if the ISM pressure rises over Pmax- 
Although this kind of mechanism can act only over short 
distances (compared with the Galaxy size) it can propa- 
gate over large scales, and the idea of self-propagated star 
formation has been used to study the formation of spatial 
patterns in galactic disks (Mueller & Arnett 1976; Gerola 
& Seiden 1978; Seiden & Gerola 1979; Schulman & Sei- 
den 1990; Jungwiert & Palous 1994). However, the phase 
transition warm gas cold clouds can be also obtained 
by decreasing Pmax under P (assumed constant) if the 
local flux of decaying neutrinos increases due, for exam- 
ple, to an increase in the supernova explosion rate. This 
triggering mechanism depends on the propagation of neu- 
trinos, and therefore can act over large distances in short 
time intervals. On the other hand, star formation can also 
inhibit in different ways the warm gas condensation self- 
regulating the star formation process inclusively over large 
distances (Cox 1983; Franco & Shore 1984; Struck-MarceU 
& Scab 1987; Parravano 1988, 1989). It has been shown 



that star formation inhibition (rather than stimulation) 
can also contribute to the formation and maintenance of 
spatial patterns in galaxies (Freedman & Madore 1984; 
Chappell & Scalo 1997). Stimulation and inhibition mech- 
anisms of star formation must be acting simultaneously in 
the Galaxy, but it is not clear yet what spatial scales are 
important for each one. The effect of non-local star for- 
mation stimulation on the formation of spiral patterns in 
galaxies should be analysed in future models. 

4. Conclusions 

The thermochemical equilibrium of the ISM, including de- 
cay of neutrinos into an ionizing photon flux F^, was cal- 
culated. The range < -Fiy < 10^ cm"'^ s^^ always shows 
two regions of stability (a warm and a cold phase) that 
can coexist in equilibrium if the ISM pressure is below a 
threshold value (Pmax)- High F^, values (~ 3 x lO'* cm^^ 
s^^) estimated by Sciama (1997a) to produce rie — 0.05 
cm~'^ at r ~ 6000 K, only can be consistent with observed 
ISM pressures if more efficient processes are heating the 
gas. It also was showed that a neutrino flux increase (due, 
for example, to an increase in the supernova explosion 
rate) may stimulate the condensation of cold gas (and 
probably the star formation) by decreasing Pmax under 
the ISM pressure value. 
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Figure captions 

Fig. la. Thermal pressure as a function of the total hydrogen 
density in equilibrium for Go = 1 and for = (solid 
Hue), = 10^ (dashed line) and = 10^ cm'^ s'^ 
(dot-dashed line). 

Fig. lb. The electron fraction as a function of the total hydro- 
gen density in equilibrium for Go = 1 and for F^ — 
(solid hne), F^ = 10^ (dashed line) and F^ = 10^ cm-^ 
(dot-dashed line). 
Fig. 2. The maximum pressure Pmax for the coexistence of 
warm and cold gas as a function of Fjj for Gq — 1 
(solid line), Go = 10 (dashed line) and Go = 20 (dot- 
dashed line). 

Fig. 3. The electron density as a function F^ for T = 6000 K 
and for Go = 1 (solid line), Go = 10 (dashed line) and 
Go = 20 (dot-dashed line). 



